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Abstract 

The dominant production mechanism for heavy quark-antiquark bound states in 
very high energy processes is fragmentation, the splitting of a high energy parton into a 
quarkonium state and other partons. We show that the fragmentation functions D(z, y) 
describing these processes can be calculated using perturbative QCD. We calculate the 
fragmentation functions for a gluon to split into S-wave quarkonium states to leading 
order in the QCD coupling constant. The leading logarithms of y/rriQ, where [i is the 
factorization scale and mg is the heavy quark mass, are summed up using Altarelli- 
Parisi evolution equations. 



Quantitative evidence for quantum chromodynamics (QCD) as the fundamental field 
theory describing the strong interactions has come primarily from high energy processes in- 
volving leptons and the electroweak gauge bosons. Such processes are simpler than most 
purely hadronic processes, because leptons and electroweak gauge bosons do not have strong 
interactions. The next simplest particles as far as the strong interactions are concerned are 
heavy quarkonia, the bound states of a heavy quark and antiquark. While not pointlike, the 
lowest states in the charmonium and bottomonium systems have typical radii that are sig- 
nificantly smaller than those of hadrons containing light quarks. They have simple internal 
structure, consisting primarily of a nonrelativistic quark and antiquark only. The charmo- 
nium and bottomonium systems exhibit a rich spectrum of orbital and angular excitations. 
Thus in addition to being simple enough to be used as probes of the strong interactions, 
heavy quarkonia are also a potentially much richer source of information than leptons and 
electroweak gauge bosons. 

In most previous studies of the production of heavy quarkonia in high energy pro- 
cesses, it was implicitly assumed that they are produced by short distance mechanisms, in 
which the heavy quark and antiquark are created with transverse separations of order l/E, 
where E is the characteristic energy scale of the process. In this paper, we point out that 
the dominant mechanism at very high energies is fragmentation, the production of a high 
energy parton followed by its splitting into the quarkonium state and other partons. The QQ 
pair are created with a separation of order 1/itiq, where uiq is the mass of the heavy quark 
Q. The fragmentation mechanism is often of higher order in the QCD coupling constant 
a s than the short distance mechanism, but it is enhanced by a factor of (E/ttiq) 2 and thus 
dominates at high energies E » itlq. The fragmentation of a parton into a quarkonium 
state is described by a fragmentation function D(z,fi), where z is the longitudinal momen- 
tum fraction of the quarkonium state and /i is a factorization scale. We calculate to leading 
order in a s the fragmentation functions D(z, uiq) for gluons to split into S-wave quarkonium 
states at energy scales fi of order itiq. The fragmentation functions at larger scales /i are 
then determined by Altarelli-Parisi evolution equations which sum up the leading logarithms 
of fi/niQ. 

One of the quarkonium processes that is important in hadron collider physics is the 
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production of charmonium at large transverse momentum pr- A charmonium state can 
either be produced directly at large px or it can be produced indirectly by the decay of 
a large pt -B-meson or a higher charmonium state with large px- In previous calculations 
of the rate for direct production of charmonium at large px Q, the dominant mechanisms 
were assumed to be short distance processes, in which a collinear cc pair in a color-singlet 
S-wave state is created with transverse separation on the order of 1/pr- A typical Feynman 
diagram which contributes to the production of the 1 So charmonium state i] c at order 
is the diagram for gg — > ccg shown in Figure 1. The order-a^ radiative corrections to this 
process include the Feynman diagram for gg — > ccgg shown in Figure 2. In most regions of 
phase space, the virtual gluons in Figure 2 are off their mass shells by amounts of order pr, 
and the contribution from this diagram is suppressed relative to the diagram in Figure 1 by 
a power of the running coupling constant ol s (jpt)- But there is a part of the phase space in 
which the virtual gluon attached to the cc pair in Figure 2 is off-shell by an amount of order 
m c . The propagator of this virtual gluon enhances the cross section by a factor of p T /m 2 c . 
At large enough p?, this easily overwhelms the extra power of the coupling constant a s . The 
enhancement is due to the fact that the cc pair can be produced with transverse separation 
of order l/m c instead of 1/pr- 

A more thorough analysis of the amplitude for gg — > r] c gg reveals that the term that 
is enhanced by Pr/m 2 can be written in a factored form. The first factor is the amplitude for 
the production of a virtual gluon g* with high-p-^ but low invariant mass q 2 via the process 
gg — ► gg*. In the limit q 2 «p\, it reduces to the on-shell scattering amplitude for gg — ► gg. 
The second factor is the propagator 1/q 2 for the virtual gluon. The third and final factor 
is the amplitude for the process g* —>■ rj c g, in which an off-shell gluon fragments into an r\ c 
and a gluon. The factoring of the amplitude allows the fragmentation contribution to the 
differential cross section da nc {E) for producing an rj c with energy E » m c to be written in 
a factorized form: 

da Vc (E) w / dz da g (E/z) D g ^ Vc (z,m c ) , (1) 

J 

where da g (E) is the differential cross section for producing a real gluon of energy E. All 
of the dependence on the energy E appears in the subprocess cross section da g , while all 
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the dependence on the quark mass m c is in the fragmentation function D(z,m c ). The 

variable z is the longitudinal momentum fraction of the i] c relative to the gluon. The physical 

interpretation of (p]) is that an r] c of energy E can be produced by first producing a gluon 

of larger energy E/z which subsequently splits into an rj c carrying a fraction z of the gluon 
energy. 

The generalization of the leading order formula ([!]) to all orders in ot s is straight- 
forward. At higher orders in a s , the gluon that splits into a quarkonium state O can it- 
self arise from the splitting of a higher energy parton into a collinear gluon. This split- 
ting process gives rise to logarithms of E/tuq. In order to maintain the factorization 
of the dependences on E and mq, it is necessary to introduce a factorization scale /i: 
\og(E/rriQ) = log(E/fi) + log(/x/mg). To all orders in a s , the fragmentation contribution 
to the differential cross section for producing a quarkonium state O with energy E can be 
written in the factorized form 

da Q (E) = V / dz ddi(E/z,fi) Di^ (z,fi) , (2) 
. Jo 

where the sum is over all parton types i. The scale [i is arbitrary, but large logarithms of 
Ej \x in the parton cross section d&i can be avoided by choosing \x on the order of E. Large 
logarithms of n/rriQ then necessarily appear in the fragmentation functions D^q(z, fi), but 
they can be summed up by solving the evolution equations 0] 

d /"l dy 

li-Q-Di^ (z,ii) = J2J — p i^A z /v^) D j^o(y,v) , (3) 
a* j z y 

where Pi^j(x, fx) is the Altarelli-Parisi function for the splitting of the parton of type i into a 
parton of type j with longitudinal momentum fraction x. For many applications, calculations 
to leading order in a s require only the g — ► g splitting function, which is 

PM - 5afei(lz- + _^_ + ^_„ + »4p*«i-„) , (4, 

where n/ is the number of light quark flavors. The boundary condition on this evolution 
equation is the fragmentation function Di^o(z,mQ) at the scale mq. It can be calculated 
perturbatively as a series in q s (wq). 
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We proceed to calculate the fragmentation function D g ^ Vc (z,m c ) for a gluon to split 
into the 1 S'o charmonium state r\ c to leading order in a s {m c ). A process (such as gg — > gg) 
that produces a real gluon of 4-momentum q has a matrix element of the form A4 a e a (q), 
where e a (q) is the polarization 4- vector of the on-shell (q 2 = 0) gluon. In the corresponding 
fragmentation process (such as Figure 2), a virtual gluon is produced with large energy 
go >> fn c but small invariant mass s = q 2 of order m 2 , and it subsequently fragments into 
an Tj c and a real gluon. The fragmentation probability Jq dzD(z, m c ) is the ratio of the rates 
for these two processes. In Feynman gauge, the fragmentation term in the matrix element 
for the rj c production has the form Ai a (—ig al3 /q 2 )Ap, where A4 a is the matrix element for 
the production of the virtual gluon and Ap is the amplitude for g* — > r] c g. The fragmentation 
term is distinguished from the short distance terms in the matrix element by the presence of 
the small denominator q 2 of order m 2 . The amplitude Ap can be written down using standard 
Feynman rules for quarkonium processes |§. Multiplying A a by its complex conjugate and 
summing over final colors and spins, we get 



where p is the 4-momentum of the r) c . Terms proportional to q a or qp are gauge artifacts and 
can be dropped. In an appropriate axial gauge, q a and qp are of order ml/q when contracted 
with the numerator of the propagator of the virtual gluon. In covariant gauges, the q a and 
qp terms are not suppressed but are cancelled by other diagrams. In the p a pp term, we can 
set p = zq + p± up to corrections of order m 2 /q Q , where z is the longitudinal momentum 
fraction and p± is the transverse part of the 4- vector p. In a frame where q = (go, 0, 0, q^), 
z — (jPo +P3)/(<?o + Q3) and p± = (0,pi,p2,0). After averaging over the directions of the 
transverse momentum, Pj_Pj_ can be replaced by —g a @p 2 /2, up to terms that are suppressed 
in axial gauge. The terms in ([5]) that contribute to fragmentation then reduce to 



J2A a A, = ( 8 -4m^ ( (3 - 4mg)2 " 2(l-*)(--4m^)H^) ■ (6) 
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We have used the conservation of the go — 93 component of the 4-momentum in the form 
s — (p± + 4m^)/z + pl/(l — z). At this point, it is easy to calculate the rate for production 
of r) c g in the limit q$ >> q 2 ~ m 2 and divide it by the rate for production of an on-shell 
gluon. The resulting fragmentation probability is 

1 a 9 n (A a " \ R ^ r f 1 * 9 s 2 + lQrni-2z(s + 4m 2 c )s + 2z 2 s 2 
o dzU ^ Z > ~ 3tt 2m c hml hml/s s 2 (s-4m 2 ) 2 ' 10 

where -R(O) is the nonrelativistic radial wavefunction at the origin for the S-wave bound 
state. We have increased the upper endpoint of the integration over s to infinity, because the 
resulting error is of order rn^/qQ, which we have been consistently neglecting. Interchanging 
orders of integration, we can read off the fragmentation function: 



1 „. ^ ^2 1^(0)1 
P 
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D n (z,2m c ) = -a s (2m c ) 2 ^f I 3z - 2z 2 + 2(1 - z) log(l - z) I . (8) 



We have set the scale in the fragmentation function and in the running coupling constant 
to = 2m c , which is the minimum value of the invariant mass y/s of the fragmenting 
gluon. In the denominator, we have set 2m c = M Vc , which takes into account the correct 
phase space limitations and is accurate up to relativistic corrections. The value of the 
S-state wavefunction at the origin -R(O) is determined from the ip electronic width to be 
|i?(0)| 2 = (0.8 GeV) 3 . We use the value a s {2m c ) = 0.26 for the strong coupling constant. 

Given the initial fragmentation function (|8]), the fragmentation function is determined 
at larger values of \x by solving the evolution equation (^) with (|j) as a boundary condition. 
The ^-dependence of D g _> Vc (z, /i) at the energy scales \i = 2m c and \i = 20m c is illustrated 
in Figure 3. The evolution causes the fragmentation function to decrease at large z and to 
diverge at z = 0. A physical cross section like (0) will still be well-behaved, because phase 
space limitations will place an upper bound on the parton energy E/z which translates into 
a lower bound on z. It is evident from Figure 3 that taking into account the evolution of 
the fragmentation function can significantly increase the rate for the production process, 
particularly at small values of z. 

The fragmentation function for a gluon into J /if} can be calculated to leading order 
in a s from the Feynman diagrams for g* — > ipgg. The square of the amplitude J2 -A.a-A.% for 



this process can be extracted from a calculation of the matrix element for e + e~ — > ipgg M. 
The calculation of the fragmentation function is rather involved and we present only the 
final result: 

,|i?(0^ 2 [» ^ /-(i-H-J/a ^ 1 



D g ^(z,2m c ) = ——^a s (2m c ) / dr / dy - — - 

IAAtt 2 Mi Jo J(r+z 2 )/2z (1 - y) 2 {y - r) 2 (y 2 - r) 2 



iff, s , s l + r-2y y - r + ^y 2 - r\ 

2^ z \Mr,y) + 9i[r,y)- f : m log ^ • 9 

t^o V 2 iy _r ) VV r y-r-^y 2 -rj 

where the integration variables are r = 4m 2 /s and y — p ■ q/s. The functions fi and gi are 

f (r,y) = r 2 (l+r)(3 + 12r + 13r 2 ) - 16r 2 (l + r)(l + 3r)y 

- 2r(3 - 9r - 21r 2 + 7r 3 )y 2 + 8r(4 + 3r + 3r 2 )y 3 - 4r(9 - 3r - 4r 2 )y 4 

- 16(1 + 3r + 3r 2 )y 5 + 8(6 + 7r)y 6 - 32y 7 , (10) 
f^r, y) = -2r(l + 5r + 19r 2 + 7r 3 )y + 96r 2 (l + r)y 2 + 8(1 - 5r - 22r 2 - 2r 3 )?/ 3 

+ 16r(7 + 3r)y 4 - 8(5 + 7r)y 5 + 32y 6 , (11) 

/ 2 (r,y) = r(l + 5r + 19r 2 + 7r 3 ) - 48r 2 (l+r)y - 4(1 - 5r - 22r 2 - 2r 3 )y 2 

- 8r(7 + 3r)?/ 3 + 4(5 + 7r)y A - 16?/ 5 , (12) 
g (r,y) = r 3 (l - r)(3 + 24r + 13r 2 ) - 4r 3 (7 - 3r - 12r 2 )y - 2r 3 (17 + 22r - 7r 2 )y 2 

+ 4r 2 (13 + 5r-6r 2 )?/ 3 - 8r(l + 2r + 5r 2 + 2r 3 )y 4 - 8r(3 - llr - Qr 2 )y 5 

+ 8(1 - 2r - 5r 2 )y 6 , (13) 

gi (r,y) = -2r 2 (l+r)(l-r)(l + 7r)y + 8r 2 (l + 3r)(l - 4r)y 2 

+ 4r(l + lOr + 57r 2 + 4r 3 )y 3 - 8r(l + 29r + 6r 2 )y 4 - 8(1 - 8r - 5r 2 )y 5 , (14) 

g 2 (r,y) = r 2 (l + r)(l - r)(l + 7r) - 4r 2 (l + 3r)(l - Ar)y 

- 2r(l + lOr + 57r 2 + 4r 3 )y 2 + 4r(l + 29r + 6r 2 )y 3 + 4(1 - 8r - 5r 2 )y 4 . (15) 

The integrals over r and y in (|9|) must be evaluated numerically to obtain the fragmentation 
function at the energy scale fi = 2m c . At larger values of /i, it is found by solving the 
evolution equation (|3|) with (g) as a boundary condition. The ^-dependence of D g ^(z,fi) 
at the energy scales fi = 2m c and fi = 20m c is illustrated in Figure 4. The evolution causes 
the fragmentation function to decrease at large z and to diverge at z = 0. 



An order of magnitude estimate of the gluon fragmentation contribution to quarko- 
nium production in any high energy process can be obtained by multiplying the cross 
section for producing gluons of energy E > 2rriQ by the initial fragmentation probability 

Jq 1 dzD(z 1 2mQ). For the r\ c and ip, these probabilities are 4.6 ■ 1CT 5 and 2.8 ■ 10 -6 . Thus, 
we can expect the asymptotic production rate of ip to be more than an order of magnitude 
smaller than r] c . The initial fragmentation function for the splitting of a gluon into the 3 Si 
bottomonium state T is also given by @, except that the mass M$ is replaced by My, the 
scale 2m c is replaced by 2rrib, and -R(O) is the appropriate wavefunction at the origin. The 
initial fragmentation probability is about 5.3 ■ 10 -7 , which is smaller than for ip by about a 
factor of 5. 

In hadron colliders, short distance processes dominate the direct production of char- 
monium at small p?, because fragmentation processes are suppressed by powers of a s (m c ). 
At sufficiently large px, fragmentation processes must dominate, because the short distance 
processes are suppressed by powers of m 2 c /p T . We can make a quantitative estimate of the pt 
at which the crossover occurs by comparing the differential cross sections for short distance 
process with the differential cross section for the gluon scattering process gg — > gg multiplied 
by the appropriate fragmentation probability. For simplicity, we consider 90° scattering in 
the gg center of mass frame. In terms of pt, we have s = Ap T - In the limit pt » m c , 
the differential cross section for gg — > gr\ c is [[1] da/dt = 81nal\R(0)\ 2 / (256M Vc p T ) . The dif- 
ferential cross section for gg — > gg is da/dt = 2437ra 2 /(128j4). To allow for fragmentation 
of either of the two outgoing gluons, we multiply by twice the initial fagmentation proba- 
bility a 2 s \R(0)\ 2 /(9nMl). The resulting differential cross section exceeds that for the short 

distance process at px ~ ^/37r/4a s M r?c « 3M Vc . The differential cross section for gg — > gip 
is [Q] da/dt = 57ra^|i?(0)| 2 M ? /,/(128p| n ) and for this case we estimate the crossover point to 
be pt ~ ^Jl.05/a s M^, « 2M^. Thus fragmentation should dominate over short distance 
production at values of pt that are being measured in present collider experiments. 

The importance of fragmentation for charmonium production in high energy processes 
can also be seen from the surprisingly large rate || for Z° — > ipcc, which is two orders of 
magnitude larger than that of Z° — > ipgg. The explanation for this is that Z° — > ipgg is 
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a short distance process, while Z° — > ipcc includes a fragmentation contribution enhanced 
by (Mz/M^) 2 . The enhanced contribution arises from the decay Z° — > cc, followed by the 
splitting c — > ipc or c — > ^c. With this insight, the lengthy calculation presented in Ref. 
|| can be reduced to a simple calculation of the fragmentation function D c ^(z, M z ). This 
calculation will be presented elsewhere IJ. 

The probability for a virtual gluon to decay into a ip was calculated by Hagiwara, 
Martin and Stirling and used to study J/ip production from gluon jets at the LEP 
collider. They did not calculate the fragmentation function D g ^(z, fi) and thus were unable 
to sum up large logarithms of Mz/m c . Furthermore, their expression for the fragmentation 
probability is missing a factor of l/(167r 2 ), which explains the surprisingly large rate that 
they found for this -^-production mechanism. 

A complete calculation of the fragmentation contribution to ip production in high 
energy processes must include the production of the P-wave charmonium states Xcj, followed 
by their radiative decays into ip. In calculating the fragmentation functions for gluons to split 
into P-wave charmonium states, there are two distinct contributions that must be included 
at leading order in a s . The P-wave state can arise either from the production of a collinear 
cc pair in a color-singlet P-wave state, or from the production of a collinear cc pair in a color- 
octet S-wave state ||. Calculations of the P-wave fragmentation functions will be presented 
elsewhere ||. 

We have shown in this letter that the dominant production mechanism for quarkonium 
in high energy processes is fragmentation, the production of a high energy parton followed by 
its splitting into the quarkonium state. We calculated the fragmentation functions D(z,fi) 
for gluons to split into S-wave quarkonium states to leading order in a s . The fragmentation 
functions satisfy Altarelli-Parisi evolution equations which can be used to sum up large 
logarithms of n/rriQ. Most previous calculations of quarkonium production have considered 
only short-distance production mechanisms, because the fragmentation mechanism is often 
of higher order in a s . At high energies E, the fragmentation mechanism dominates because 
it is enhanced by a factor of (E/itiq) 2 . In the case of charmonium production at large pt in 
hadron colliders, we estimated the transverse momentum at which fragmentation begins to 
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dominate to be less than 10 GeV. All previous calculations of quarkonium production from 
this and other high energy processes must therefore be reexamined, taking into account the 
possibility of production by fragmentation. 
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Figure Captions 

1. A Feynman diagram for gg — > ccg that contributes to r\ c production at order c? s . 

2. A Feynman diagram for gg — > ccgg that contributes to i] c production at order a A s . 

3. The fragmentation function D g ^ Vc (z, /i) as a function of z for /i = 2m c (solid line) and 
\x = 20m c (dotted line). 

4. The fragmentation function D g ^(z,ii) as a function of z for /i — 2m c (solid line) and 
ji = 20m c (dotted line). 
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